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We present here a set of lecture notes on quantum systems with time-dependent bound¬ 
aries. In particular, we analyze the dynamics of a non-relativistic particle in a bounded 
domain of physical space, when the boundaries are moving or changing. In all cases, 
unitarity is preserved and the change of boundaries does not introduce any decoherence 
in the system. 
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1. Introduction 

We present here the notes of three lectures given by one of us at the International 
Workshop on Mathematical Structures in Quantum Physics, held in February 2014 
in Bangalore at the Centre for High Energy Physics, Indian Institute of Science. 
The course considers some aspects of quantum systems with time-dependent bound¬ 
aries, a very active area both from the mathematical point of view, see for instance 
the works of Yajima m. Dell’Antonio et al [3] and Posilicano et al BE], and from 
a physical perspective. Notable applications arise in different fields ranging from 
atoms in cavities m to ions and atoms in magnetic traps [5], to superconducting 
quantum interference devices (SQUID) to the dynamical Casimir effect in mi¬ 
crowave cavities |10| . Moreover, the role and importance of boundary conditions at 
a fundamental level have been stressed in an interesting series of articles, see 
and references therein, where varying boundary conditions are viewed as a model 
of space-time topology change. 

These notes will be organized as follows. In the first lecture, in Section we 
will study the problem of a quantum particle moving in a one-dimensional box. 
In particular we will introduce and describe the problem in a pedagogical way 
focussing on the mathematical aspects of the physical model. We will try to give 
all the mathematical tools needed to a deep understanding of the physics of the 
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problem. Moreover, at the end of the section, we will follow Berry m and show 
how, quite surprisingly, even an easy spectral problem like this can lead to intricate 
and rich dynamics. In the following lectures, we will look at the dynamical evolution 
with time-dependent boundaries. In particular, in the second lecture, in Section]^ 
we will discuss the case in which the walls of the bounding box are moving mi, 
while in the last lecture, in Section we will look at the problem of the change 
of boundary conditions with time m- Finally, in Section we will draw some 
conclusions. 

Interspersed between the lectures we offer a selection of exercises ranging in 
difficulty. They are useful to focus the attention on specific points and to test the 
learning of the reader, therefore we kindly recommend not to skip them. 

2. Lecture 1: A particle in a one-dimensional box 

One of the very first problems involving boundary conditions in quantum mechanics 
is the study of a non-relativistic particle in a one-dimensional box. Indeed, it is one 
of the exercises sometimes professors give to students just to test their ability. 
The problem consists in a massive particle moving freely in an infinitely deep well 
located, say, at / = [a,b] with b > a, from which it cannot escape (Fig. [^. 



Fig. 1. A particle moving in a one-dimensional box. 


We will see that even if the problem is apparently fairly simple, it leads to 
very interesting results that stress the difference between classical and quantum 
mechanics. In fact, the very presence of boundaries in general tends to enhance 
the quantum aspects of the system. The reason is that in quantum mechanics 
the behavior at the boundary is constrained by the very structure of the theory 
and in particular by unitarity. This is at variance with classical mechanics, where 
the behavior at the boundary does not derive from fundamental principles, and is 
usually treated phenomenologically. 

From the mathematical point of view the system is described by a Hamiltonian 
operator associated to the kinetic energy of the particle, namely the Laplacian, 

r= —= (1) 

2m 2m dx"^ ’ 

acting on square integrable functions (remember that \ip{x)\'^ is a probability den- 
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sity!) with square integrable second (distribution) derivative on the interval / of 
the box, = {^|} € \ G the mathematical literature 

is known as the second Sobolev space; it is nothing but the maximal domain of 
definition of the kinetic energy operator T. 

Equation Q describes the action of T in the bulk of the box. In order to 
get a well-defined dynamics, Eq. Q should be equipped with suitable boundary 
conditions, which specify the behavior of the particle at the walls. Namely, the 
wave functions tp to which we are allowed to apply T must belong to the domain 

D{T) = {^ £ 'H^(/),with suitable boundary conditions}. (2) 

The question that here arises naturally is how to choose these boundary conditions. 
This marks the difference between classical and quantum mechanics: in quantum 
mechanics the behavior at the boundary, encoded in the domain D{T), derives 
from basic principles. Indeed, in quantum mechanics the conservation of probability 
leads to the necessity for the evolution to be represented by unitary operators. But 
unitarity of the evolution is equivalent to self-adjointness of its generator, as stated 
in Stone’s theorem m- 

Theorem 1. (Stone) Let {Ut)teR be a strongly continuous one-parameter group of 
unitaries on a Hilbert space %. Then there is a unique self-adjoint operator H on 
TL, called the generator of the group, such that Ut = for all t. 

According to Stone’s theorem the operator T, that generates the dynamics of 
the particle, Ut = must be self-adjoint, i.e. 

T = T^ and D{T) = D{T^). (3) 

In order to find the domain D{T) let us consider (j),ip € TU{[a,b\)-, then a double 
integration by part gives (prove it!): 

-{’(PW) = -[ ip{x)(p”{x)dx =-{-(("lei))-\- (4) 

J a 

A{ip, (p) = f}'(h)(t>{b) - f}'{a)(t>{a) - il}{b)cj)fb) -k 'ip{a)cj)fa). (5) 

The self-adjointness requirement amounts to finding a maximal subspace D{T) of 
TU{[a,b]) on which the boundary form A identically vanishes, namely A{ip,(j)) = 0 
for all V' € OiT). 

The more common boundary conditions you can encounter in the literature are: 

(1) Dirichlet boundary conditions: 

ip{a) = ip{b) = 0; (6) 

(2) Neumann boundary conditions, involving the derivative of the function: 

V’'(a) = i’'{b) = 0. (7) 

(Of course, the above conditions must be satisfied by both tp and (p in the 
boundary form A{tp, (p)\) 
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Exercise 1. Check that the Dirichlet and the Neumann boundary conditions make 
A identically vanish, and thus, modulo a check of maximality, yield a good domain 
D{T) of self-adjointness for the kinetic energy operator T. 

There is a one-parameter family of boundary conditions that includes Dirichlet and 
Neumann as limit cases: Robin’s boundary conditions, that connect the function 
with its derivative at the boundary: 

1 Q! 1 q; 

V''(a) =--tan-?/'(a), tp'(b) = - tan - tjjib), (8) 

Lq Z Lq Z 

where a £ (—7r,7r], and Iq is a reference length. From Q we recover Dirichlet and 
Neumann in the limits a —>■ tt and a —>■ 0, respectively. 

Notice that the choice of a domain of self-adjointness D{T) is not just a mathe¬ 
matical nuisance. Different domains give rise to different physics! Indeed, the choice 
of boundary conditions strongly depends on the physical behavior of the walls, as 
we will now show. 

Let us concentrate on one wall, e.g. the left one at a = 0, and look at the 
scattering process of a plane wave arriving from the left with momentum —hk < 0. 
After reflection by the wall a plane wave with momentum hk will arise. Therefore, 



Fig. 2. Reflection at a wall 


the wave function will have the form 

V^(a:)~e-‘'=^-bre*'=“, (9) 

where r is a reflection coefficient. Since we are dealing with an impenetrable box, 
and we assume that the particle does not stick at the wall, we are forced to require 
that |rp = 1, so that the only freedom left is in choosing a phase, namely r = e'^. 
The Dirichlet boundary condition, ^(0) = 0, corresponds to a phase /3 = tt for all 
k, while Neumann corresponds to the choice /3 = 0 for all k. Notice that in general 
the phase f3{k) has a nontrivial dependence on the wavenumber k. The following 
exercise shows that for a generic a in (|^ the behavior of the walls depends on the 
momentum hk of the impinging particle through a nontrivial phase f3{k). 

Exercise 2. Prove that Robin’s boundary condition ~ tan ^ ^/>(0) corre- 
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spends to a momentum dependence of the phase shift /3{k) given by 


tan 


m 


1 a 

—— tan —. 

k Iq 2 


( 10 ) 


Recover, as a limit, Dirichlet and Neumann, and show that they are the only phases 
independent of k. 


The general theory of self-adjoint extensions and boundary conditions 
asserts that all possible boundary conditions that make T self-adjoint can be 
parametrized by 2 x 2 unitary matrices U in the following way: 


where 


i(7 + C/)T' = (J 



; = 


h 


u eu{2), 


( 11 ) 

( 12 ) 


and I is the identity matrix. 


Exercise 3. Check that if ■0 and (j) satisfy for the same unitary C7, then 

the boundary form ([^ vanishes, A('0, 0) = 0. 

Robin’s boundary conditions (|^ correspond to (prove it!) 


U = e-'“J = 



(13) 


and in particular the Dirichlet boundary conditions correspond to U = —I, while 
the Neumann ones are given by 17 = I. More generally, by taking 


U = 



(14) 


with ai,a2 G (—7r,7r] we can describe different boundary conditions at the two 
endpoints, e.g. Dirichlet/Neumann (ai = tt, a 2 = 0), and include all cases con¬ 
sidered above, when ai = 02 - But, of course the unitaries (141 do not exhaust all 
allowed boundary conditions, since they form only a two-dimensional submanifold 
m X 17(1) of the total four-dimensional manifold 17(2). Thus, where do all missing 
boundary conditions come from? 

Up to now we have considered only local boundary conditions, somewhat de¬ 
ceived by the physics depicted in Fig. the two walls did not talk each other. 
Indeed, the unitary matrices in (14) are diagonal and do not mix the boundary 
values at a with those at b; we missed all nondiagonal unitaries, describing nonlocal 
boundary conditions. Mathematics tells us that unitarity allows also for them. So 
what physics do they describe, if any? 

In fact, nondiagonal unitaries are describing a physical situation which is dif¬ 
ferent from that of a box with two walls far apart. In order to interact, the two 
ends of the interval should come close, as in Fig. so that the interval should be 
bent and the two walls should become the two sides of a junction. In other words. 
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Fig. 3. Change of topology: from an interval to a circle 


by changing from a diagonal 17 to a nondiagonal U we are assisting at a change of 
topology: from an interval to a circle. 

Therefore, the geometry that is able to support all possible boundary conditions 
in hl{2) is that of a ring with a junction. If the junction is impermeable, i.e. there is 
total reflection at the walls, we are back to the interval, otherwise there is nonzero 
transmission across the junction, from one wall to the other. An interesting example 
is given by the matrix 


U = 


/O e-‘“\ 

Ve‘“ 0 ) 


= cos aax + sin a Uy , 


(15) 


which describes pseudo-periodic boundary conditions (prove it!): 
ii){b) = e‘“ V'(fl)j V'^(^) = 


(16) 


By passing through the junction, the wave function acquires a phase a. If a = 0 
we get the famous periodic boundary conditions, and the geometry is that of a 
circle; if a = tt we get antiperiodic boundary conditions. The phase a encodes the 
properties of the junction, for example the material it is made of, or its width. 

In general, if the unitary U has no —1 eigenvalues, the wave function ij) can 
assume any value at the endpoints of the interval. Only the boundary values of its 
derivative are constrained in some way. On the other hand, one eigenvalue equal to 
— 1 corresponds to one constraint on the values of 'ip at the ends, as for example in 
the first equation in (16). Finally, two —1 eigenvalues, i.e. U = —I, correspond to 
two constraints on the boundary values of the wave function, i.e. Dirichlet at both 
ends. 


Exercise 4. Prove that the unitary (15) has always an eigenvalue equal to —1. Find 
its corresponding eigenvector ^ and show that the first equation in (161 is nothing 
but an orthogonality condition (^I'h) = 0. (In Lecture 3 we come back to the 
geometrical meaning of this condition.) 


2.1. Hard walls 

Let us concentrate now on the textbook case of a particle confined in an infinitely 
deep well. Fig. The appropriate boundary conditions are Dirichlet’s: '4’{a) = 
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tjj{b) = 0. The dynamics of the particle is described by the Schrodinger equation 

.^dip{x,t) _ b? 

' at 2m dx^ ' ^ ^ 

A separation of variables, ip{x,t) = u{x) exp{—iEt/h), reduces the problem to the 
solution of the spatial part of the differential equation, that means to find eigen¬ 
vectors and eigenvalues of the operator T: 

— - — u" {x) = Eu{x). (18) 

2m 

The general solution is 


u{x)=cie+''’^^+ C2e-'^^^, (19) 

with k = yj2mE/h (in principle k can be imaginary, but see below), and Ci and C 2 
are arbitrary constants that can be fixed (up to a common phase) by imposing the 
Dirichlet boundary conditions, 


u{a) = Cl e*'=“ -b C 2 e-‘''“ = 0 (20) 

m( 6) = Cl + C2 = 0, (21) 

and normalization 


{u\u) = 1. 


( 22 ) 


Exercise 5. Prove that the normalized eigenfunctions of T, with the Dirichlet 
boundary conditions, are 


Unix) = 



(23) 


where I = b — a, and that the eigenvalues, giving the permitted energy levels are 


for n = 1, 2,-See Fig. 


En 


b? n'^Tr'^ 
2m P ’ 


(24) 


Remark 1. Eq. (24) tells us two important things on the bound states of the 
particle: 


(1) the energies are quantized; 

(2) the energy is always strictly positive. 


The strictly positivity of the energy is a property of the Dirichlet boundary 
conditions. Indeed, with the Neumann ones the energy of the ground state is 0, 
and, more surprisingly, it is even negative with Robin’s boundary conditions! 







8 S. Di Martino, P. Facchi 



Fig. 4. The eigenfunctions (solid lines) of the operator T in the box at different levels of energy 
(dashed lines). 


Exercise 6. Prove that the ground state of T, with the Neumann boundary condi¬ 
tions is 


and has zero energy, Eq 
boundary conditions Q. 


Vo{x) = 



(25) 


0. Then, look at the eigenvalue problem with Robin’s 


2.2. Fractals in a box 


The textbook exercise of the quantum particle in a box inevitably ends with the 


evaluation of the eigenvalues (241 and the eigenfunctions (23). The result is so simple 


and intelligible that we all felt a profound satisfaction when we derived it in our first 
course of quantum mechanics. The simplicity of the spectrum is deceptive and leads 
us to think that we fully understand the physical problem. In particular, we are 
convinced that the dynamics, which is the solution to the Schrodinger equation 0, 
must surely be as much simple. In fact, this belief is false, as showed by Berry m- 
the dynamics is instead very intricate. 


Let us assume at time t = 0 that 'ijj{x,0) = vo{x), with vq given by (25). This 


is the simplest conceivable initial condition, corresponding to a flat probability in 
the box [a, 6], with I = b — a. We are interested in the time evolution of this initial 


wave function. Its L^-expansion in terms of the eigenfunctions (23) of T reads 


OO 


Xo(.x) = > CnUn{x), 


( 26 ) 
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Fig. 5. Graphs of |'i9(Gt)P vs ^ at different rational times r along the Fibonacci sequence tending 
to the golden mean, ^ = (1 + ■v/5)/2. See the emergence of a fractal structure. 


where c„ = {un\vo). 

Exercise 7. Show that 

= (27) 

and, in particular, C 2 „ = 0, for all n = 1, 2,.... 


In the same way the quantum evolution, described by the action of the unitary 
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operator U{t) = e 

ip{x,t) = ^ c„ Un{x). 


, can be written as an L^-convergent series 
+ 00 


n—1 


We can now use the explicit expressions (23l-(24) and obtain 

+ 00 


ip{x,t) = Jj Cn sin (^ (x - g) j exp 


ih n^TT^t 
2m P 


(28) 


(29) 


In terms of the dimensionless variables ^ = I ^{x — {a + b)/2) S [—1/2,1/2], and 


T = 2tt th/mr S K, it reads 
+ 00 

' C 2 k+i sin 


= Vy E 


k^O 


2^(^+2 




exp 


—inr 1^+2 


(30) 


By writing the sine as the sum of exponentials and by making use of the expres¬ 
sion (27), we finally get 


+ O0 


d(e,r)= 

n——oo 

(notice that now the sum runs over all n G Z). 


Ce 


1 1 
2’ 2 


(31) 


Exercise 8. Derive Equation (31) and show that 

1 (- 1 )" 


dn — 


X^/l ' 


1 ’ 
2 


n € Z. 


(32) 


If we take a closer look at the expression (31) we notice that it is a Fourier 
series with quadratic phases. This series is the boundary value of a Jacobi theta 
function [18) . which is defined in the lower complex half plane of r, and it has a 
very rich structure investigated at length by mathematicians. For a full immersion 
in its deep arithmetic properties see the charming “Tata Lectures on Theta” by 
Mumford [19j . A simple property is its quasi-periodicity in the rescaled time r 
(check it!): 

d(C,T + l) = e-'"/4d(^,r). (33) 

Thus at integer times t the wave function comes back (up to a phase) to its initial 
flat form (25): these are the quantum revivals. More generally, at rational values of 
T the graph of |'d(^, r)p is piecewise constant and there is a partial reconstruction 
of the initial wave function [20], see Fig. On the other hand, at irrational times, 
the wave function is a fractal, with Hausdorff dimension Dh = 3/2, as shown in 
the last panel of Fig. In fact, 'd(^, r) can be proved to be a fractal function in 
space and time, and to form a beautifully intricate quantum carpet, with different 
Hausdorff dimensions along different space-time directions |21j . 
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3. Lecture 2: Moving walls! 

In this lecture we will try to answer the following question: What happens if the 
walls of our box start moving? Which equation will describe the quantum dynamics 
of the bouncing particle? 

The classical version of this problem was introduced by Fermi [22] in 1949, and 
then investigated by Ulam |23j . It is convenient to parametrize the conhning interval 
as 

1 j 

—, d “t“ — 

2 2 


Il,d — 


d- 


(34) 


so that / > 0 is the width of the box and d S M its center. We will suppose that 
I and d are regular functions of time, t I{t) and t i—)■ d(t), with l{t) > Iq >0 
so that the interval never shrinks to a point. From now on we will often omit the 
dependence on t. 

As we did in the case of still walls, we analyze the dynamics described by the 
Schrodinger equation 




where the domain of is 


A.d = |'0 e “ 0 ='ip[d- 


= 0 


(35) 


(36) 


(Dirichlet’s boundary conditions). Notice that this domain depends on time, so that 
at different times we work on different spaces. This means that the time derivative, 

-ipit) 


d , , , tbU - 

= lim -- 

dt £->-0 


(37) 


involves the sum of vectors belonging to different Hilbert spaces, since in general 
Di{t),d{t) 7^ Therefore, we need to take more care in the formulation 

of the problem and in the interpretation of Eq. (351. 


The correct formulation of the problem can be accomplished by embedding the 
space of square integrable functions on the interval, in the larger Hilbert 

space L^(IR) on the real line: 




(38) 


where = K\X denotes the complement of the set X. Thus, every wave function 
Tp £ L^(IR) can be written as a sum i/; = y + c/), where y £ and (j) £ L^(/f^). 

Following this strategy the Hamiltonian of the system, representing the kinetic 
energy of the particle in the box, in the containing space reads: 


= I) = -£5?®“ »■ 


(39) 


equipped with the Dirichlet boundary conditions. Notice that the choice of the 
Hamiltonian to be 0 in the component L^(/f^), i.e. outside the box, is completely 
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arbitrary and immaterial, since we will be interested in the dynamics inside the 
box, for a particle with initial (and evolved) wave function with (j) = 0. Moreover, 
it is worth noticing that even if the Hamiltonian in Eq. (35) appeared to be time- 
independent, the direct sum decomposition in (391, in the case of moving walls, 
makes its time dependence clear. Indeed, the direct sum in Eq. (38) depends on time 
through the functions I and d. In other words, the restriction to the interval somehow 
concealed the time-dependence of the Hamiltonian, which becomes evident once one 
enlarges the space. Therefore, the Schrodinger equation we are dealing with has a 
time-dependent Hamiltonian on a time-dependent domain. 


3.1. Reduce and Conquer 

The attack strategy that we will follow is to find a unitarily equivalent problem 
where the Schrodinger operators act on a common fixed domain of self-adjointness. 
Let us consider the transformation 

U{l,d) : L‘^{R) ^ (40) 

that acts as 

= + (41) 

where > 0 is a reference length. 


Exercise 9. Prove that the transformation U{l,d) is unitary and maps the interval 
Ii^d onto the reference interval 


I = 


dlo,0 


lo lo 
2' 2 


(42) 


In fact, U(l,d) can be rewritten as a composition of two transformations 


where 


U{l,d)=V 



{T{d)ip){x) = ipix 


^)V(d)t=I?(-ln^)r(-d), (43) 

- d), (l?(s)'0)(a;) = e“®/^'0(e“®a;), (44) 


are two (strongly continuous) one-parameter unitary groups implementing the 
translations and dilations on L^(IR), see Eig.[^ 


Exercise 10. Prove the decomposition (Hi 


Now notice that the domain (36) is mapped onto the fixed domain D = 

U{fd)Di^d 




( 45 ) 
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1 

d 






Translation 




Fig. 6. The effect of the transformation on the box can be described as the composition of a 
translation and a dilation. 


describing the Dirichlet boundary conditions in a fixed box /, and the Hamilto¬ 
nian (39) is mapped into 


2 2 
P 


//(O = t/i/.t/'= (I) = 


2m 




(46) 


Exercise 11. Prove that the unitary transformation U{l,d) yields constant boundary 
conditions only when the latter do not mix the wave function with its derivative. 

Henceforth we will denote the wave functions in the frame with moving and 
fixed walls by '^(a:) with x G Ii^d and by with i G I, respectively. By deriving 
the relation 


= U{l{t),d{t))ip{t), (47) 

we get the Schrddinger equation in the fixed reference frame from the one in the 
frame with moving walls: 

i/i^^ = ih (U{1, d)ip + U{1, d)V') 

= (^U{I, d)Ho{l, d) + ihil{I, d)j 4’ 

= (u{l,d)Ho{l,d)U^{l,d)+ihU{l,d)U\l,d)'jUil,d)4;. (48) 
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Notice that now the Schrodinger operator contains not only the transformed Hamil¬ 
tonian H{1), but also an additional geometrical term 


K(l,d) =ihU{l,d)U^l,d). 


(49) 


Let us compute this geometric contribution step by step. First of all the action of 
U = dU{l{t),d{t))/dt on a test function ip: 


dU 


dt 




dt 


I 


I 


2-\/l(jl \ I 


I 


I / I 




h \ ^0 


I 


Then, since 


we have 




(50) 


{U\l,d)(j)){x) = iTid)V{\iil/lo)(j)){x) = ~ d)J > (51) 


(52) 


that is 


= -- [xp — i-^ \ — '-ydp. 


dt 


I 


I 


(53) 


with X and p the position and momentum operators. Thus, the geometric generator 
of the unitary transformation reads 


K{1, d) = iH = --xop — 


(54) 


where AoB = {AB + BA )/2 is the symmetrized (Jordan) product of the operators 
A and B, and the canonical commutation relation [x,p] = ih has been used. 


Exercise 12. Show that the generator of the translation group T{d) in (44) is the 
momentum operator p, while the generator of the dilation group 2J(s) is the virial 
operator a; op. As a consequence, by using the decomposition (43) show that 

U{I, d) = exp ( ^ In ( ^ ) X o p ) exp ( ) , (55) 


which yields (54). 

Summing up, we finally come to the Schrodinger equation in the reference frame 
with fixed walls 


= (H(i) + K{l.d))4, = 


( 56 ) 
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Let US now make some considerations about this equation and its solution. We 
need a definition and a theorem m- 

Let T and A be two operators on a Hilbert space 'H, such that D(T) C D{A). 
A is said to be relatively bounded with respect to T or simply T-bounded, if there 
exist two non-negative constants a and b such that: 

IIHV^II <a||i/;||+6|lri/>||, yiPeD{T). (57) 

Moreover, the infimum of the possible values of b is called the T-bound of A. 

Theorem 2. (Kato-Rellieh) Let T be self-adjoint and bounded from below. If A is 
symmetric and T-bounded with T-bound smaller than 1, then T -\- A is self-adjoint 
and bounded from below. 

Exercise 13. Show that K{l,d) is relatively bounded with respect to H{1) with 0 
relative bound. 


Applying the Kato-Rellich theorem it is easy to prove that the total Hamiltonian 
i7(Z) -I- K(l,d) with domain D{H{1) -\- K{l,d)) = D{H{1)) — D is self-adjoint. 
Therefore, the Schrodinger equation in (561 is well defined for any initial condition 
4>{0) G D, and yields a unique unitary propagator. The proof is a corollary of 
Theorem X.70 in Plj, since t i—)■ H{l{t)) -\- K(l{t),d{t)) is a one-parameter family 
of Schrodinger operators on a common domain of self-adjointness D for any pair of 
differentiable functions d{t) and l{t), with l{t) > l^j, for some /g > 0- 


Exercise 14. Prove that the energy rate equation of the particle is given by 

to 


^2 / I 


dt 


^uit)\Hm)m = -^{T^ 


2m \l{t) 




(58) 


(59) 


Exercise 15. [Rigid translation] Find the Schrodinger operator in the moving refer¬ 
ence frame in the case of a translation without dilation of the walls: 


where t) = 9j(/)(^,f). Hint: The energy rate can be computed as 

m = - (KHcficf)) + (m^), 

where K{1, d) is the operator in (|54|. Pay attention to the domains. 


I = 0 , 


d{t) = do -I- vt. 


(60) 


3.2. An Example: The Accelerating Box 

Let us consider the case in which the box is moving with a constant acceleration g, 
without dilating: 

/ = 0, d(t) = do - 


(61) 
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Since the particle is confined in a accelerating rigid box this is the problem of a 
particle in a rocket, see Fig. If we assume I = Iq we find: 

H=^+gtp, (62) 

with Dirichlet boundary conditions. We can now apply the gauge (unitary) trans- 



Fig. 7. The particle in an accelerating box. 


formation 


G{t) : L\I) ^L\l)-.c^^x = cm, (63) 

with 

xiO = iCmiO = , (64) 

in such a way that the Hamiltonian becomes 

G{t)HG^{t) = (65) 

Therefore, the Schrbdinger equation reads 

describing a particle in a constant gravitational field, in agreement with the equiv¬ 
alence principle. 

2 

Exercise 16. The eigenfunctions of the operator ^ — mgx belonging to the eigen¬ 
value E are linear combinations of the two Airy functions [18| : 


^{x) = ciAi c^Bi f y - . 

\ t mgl J \ i mgL J 

Prove that the permitted values of energy are the solutions of: 

\ 2 mgl J \2 mgl J \2 mgl J \ 2 mgl J 

and in particular that the first four are (e = ;^), see Fig. 8 


(67) 


( 68 ) 
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ei 

9.86851 

£2 

39.4787 

£3 

88.8266 

£4 

157.914 



Fig. 8. The zeros of the function (|68[l yield the permitted energy levels. 


4. Lecture 3: Changing boundary conditions 


In this lecture we will consider a somewhat different situation: we will assume 
that the walls of the box are fixed, so that / = [a, b] is given, but the way the 
walls interact with the quantum particle changes in time due to some physical 
mechanism. That means in our model that the boundary conditions change in 
time. More generally, the situation we have in mind is that depicted in Fig.[^ a ring 
with a junction whose properties are time dependent, so that all possible boundary 


conditions (111 can in principle be implemented. The evolution of the particle will 


be given by a Schrodinger equation with a time-dependent Hamiltonian: 

with 'ip{0) = '00- Here = p'^l2m on V?{I) C with boundary condi- 


(69) 


tions (11), with U = U{t). Thus the Hamiltonian depends on time through its 


boundary conditions U(t). 

This model can be implemented in a superconducting quantum interference 
device (SQUID) with a tunable junction, obtained by replacing the junction with 
an additional flux loop |25l26l27j . This can be an experimental realization of a 
continuous change among different topologies [12] . 

We will not study this problem in full generality, but instead we will consider 
the particular case of boundary conditions rapidly alternating between two values 
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U and V (for example they periodically alternate between Dirichlet and Neumann): 

V, 2kT < t < {2k + l)r, 


U{t) = 


[/, {2k + l)T<t< {2k + 2)r, 


(70) 


with A: = 0,1, 2,, and r the time period. 

In this scenario, since the Hamiltonian is constant during the period among 
the changes, the Schrodinger equation (69) can be explicitly integrated (the trick 
is revealed: the reason for the funny time-dependence of U{t) is just this!) and at 
time 2Nt we end up with 

iI}{2Nt) = (e-irTu/h^-irTv/n'j ^^.irTu/h^-irTv/h'j ^ ^ /h^-irTv/h'^ 


N times 


= e 


-irTu jh -irTv/h 


)%o. 


(71) 


We want to study the behavior of the dynamics in the limit iV —> oo, when the 
time interval between the switches t = t/N goes to zero, the number of switches 
goes to infinite, while the total time of the evolution 2Nt = 2t is kept constant. 
Therefore, we are led to the study of the limit of the product formula: 

/ \ N 

lim . (72) 

N^oo V / 


Exercise 17. Prove that, given two matrices A and B the Lie product formula holds 


(^^-iA/N^-iB/Ny 




with C = A + B. 

Hint: use the telescopic equation 

Af-l 

-E^ =Y^ D^{D - E)E 


^N-l-k 




with D = e iS/w E = e dA+B)/N^ together with the estimates 


D-E = Q-^AtN^-iB/N _ ^-i(A+B)/N ^ _ 


2iV2 


[A,B]+0 


1 

7V3 


(73) 


(74) 


(75) 


Trotter in [28] proved that Lie’s product formula (73) is still valid when applied 
to (unbounded) self-adjoint operators A and B, such that their sum C = A + B is 
still (essentially) self-adjoint on the intersection of their domains, D{C) = D{A) n 
D{B). 

Well, this seems to be our case: tTu/h and tTy/h are unbounded self-adjoint op¬ 
erators, so the limit of the Lie-Trotter product formula (72) should be exp(—i(T£/ -|- 
Tv)t/h). Now we get 


T + T - - 2 


(76) 
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on the common domain D{Tu) n D{Tv). Therefore, 

(Q—'itTu/hN^—itTv/UN \ y ^—\2tTw/fi ( 77 ) 

when iV —> oo, and we arrive at the intelligible result that rapidly alternating free 
evolutions with two boundary conditions U and V yield again a free evolution, 
possibly with different boundary conditions W. 

Unfortunately, life is not so easy: Trotter’s theorem does not hold, because 
the intersection of the the domains D(Tw) = D{Tu) n D{Tv) is too small, being 
defined by too many constraints (those of U and those of U). For example, let us 
consider Dirichlet, U = —J, and Neumann, U = /, boundary conditions. In this 
case D{T-j) n D{Ti) = {■0 | dt = d/' = 0} and in this domain the kinetic operator 
/2m is symmetric but not self-adjoint! In other words, our common sense tells us 
that we should get /2m with some boundary conditions VF, but we do not know, 
even heuristically, which ones, because the boundary conditions U and V fight each 
other and yield too many constraints! 

Math again comes to our aid: The right objects to look at are not the operators 
Tjj and Ty, but the quadratic forms associated to them, that is their expectation 
values tu{ip) = {if/Tui/!) and ty(V') = {4’\Tv4’)i measuring the kinetic energy of 
the particle in state ip. The crucial fact is that the quadratic form a{ip) = {ip\Aip) 
associated to an unbounded operator A makes sense on a domain D(a) which is 
larger than the domain of A, namely D{A) C D{a). Thus the sum of two quadratic 
forms a{ip) + b{ip) is well defined on wave functions ip G D(a) n D{h) for which the 
operator sum A -|- B is not defined, since ip ^ D{A) n D{B). 

The extension of Trotter’s theorem to quadratic forms reads as follows [29130]: 


Theorem 3. (Lapidus) Let a and b be the quadratic forms associated to A and B. 
If A and B are self-adjoint and bounded below, and D — D{a)r]D{b) is dense, then 
Trotter’s formula (73) holds with 


C = A+B, 


(78) 


the unique operator associated to the quadratic form a-\-b. The operator C is called 
the form sum of A and B. 


In our case, we will see that the domains of the quadratic forms of the kinetic 


energy tij{ip) involve only the boundary values (121 of the wave function dt and not 


the boundary values of its derivative d'h This will imply that the sum of the kinetic 
energies 


twi'f’) = 


tuilp) +tviip) 


(79) 


is always the quadratic form associated to a self-adjoint operator T^, namely 

tw{ip) = i'f’lTw'^P), i’ G D{Tw). (80) 
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For example, in the case of Dirichlet and Neumann we will see that D{t-i) = 
I = 0} and D{ti) = {i/i | no bound, conds.}, so that D{ti) = {^ | v]> = 

0} = alternating Dirichlet and Neumann give Dirichlet. 


Summing up, we get that the limit (771 holds with 


Tw = 


Tu + Tv 


(81) 


the form sum of Tu and Ty ■ 


4.1. From operators to quadratic forms 

For any '0 € one integration by parts yields 

.2 rb 




= £(IW'f-W» 


(82) 


where and 4*' are the two-dimensional vectors of boundary data defined in (12). 


Exercise 18. Prove that Eq. (82) holds. 


Notice that, at variance with our starting point, the last line of (82) involves 


only the first derivative of 0, and thus makes sense on a larger domain. Moreover, 
by using the boundary conditions (l^-(@ we will trade the boundary values of the 
derivative for the boundary values of the function in (82), obtaining the following 
expression for the quadratic form associated to Tu'. 


tuW = ^ (110 


/||2 


TuW), 


(83) 


with r(7(4') a quadratic form of the boundary vector dr. 

The form tu{'4^) is well defined for any integrable function with square integrable 
first (distribution) derivative, 77^(7) = {0 € T^(7) | p0 S i^(/)} (the hrst Sobolev 
space). In order to get the explicit expression of Fj/ and the precise domain D{tu) we 
have to distinguish among three possibilities according to the number of eigenvalues 
Ui and U 2 of U equal to —1. 


(1) If both the eigenvalues of U are different from —1, such as in the case of the 
Neumann boundary conditions, then (/ -|- U) is invertible, and the boundary 
values of the derivative can be expressed in terms of the boundary values of the 
function: 


4-' = 


i I-U 
hiFU 




(84) 
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Therefore we get 




(85) 


and D{tu) = W{I), with no constraints on the boundary values 
(2) If U has only one eigenvalue equal to —1, namely ui = —1 and U 2 ^ —1, then, 
if we call ^ the normalized eigenvector corresponding to the eigenvalue —1, and 


its orthogonal, from (11) we get 


(CI«') = o, (ri^') = - 




Iq 1 U2 




Therefore, 


(eivi/) = o, = 

io 1 + U 2 


( 86 ) 


(87) 


Exercise 19. Prove Equations (86) and (87). 


(3) Finally, if ui = U 2 = — 1, i.e. in the case of the Dirichlet boundary conditions, 
then U = —I, so that 


= 0 , 


r-7(^) = 0. 


( 88 ) 


4.2. Composition law of boundary conditions 


Now we can use the quadratic form (83) to evaluate the limit of the alternating 


dynamics 0, according to the recipe fM). By cooking our equations we will prove 
that the boundary conditions W are given by the composition law: 


w = u*v = v*u, 


(89) 


where * is a commutative and associative product on the boundary unitary opera¬ 
tors. 


The evaluation of the emergent boundary condition W in (89) requires the 


computation of the sum of the kinetic energies (79) and the evaluation of the domain 

D{tw) = D{tu) n D{tv)- (90) 

Again, we distinguish various cases according to the number of eigenvalues —1: 

(1) When both U and V have no eigenvalues equal to —1, we get 

Tw = iTu + rv)/2, (91) 

with no constraints on the wave-function boundary values 'k. By using the 


expression (85) we obtain (prove it!) 


w = u*v = 


I_lfCzE + Czv) 
2 [l+U ^ I+V I 


I+^(^ 

^ ^ 2 \I+U 


l-U 


I-V 

I+V 


(92) 
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(2) If —1 is a nondegenerate eigenvalue of U, and V has no eigenvalues equal to 
— 1, then D[tw) = D{tu), with the only constraint (Cl'h) = 0. Therefore, the 
boundary forms T;/ and Ty are nonzero and add up only on the orthogonal 
subspace, spanned by It is easy to see that 

IT = C/*y = -|e)(e|+u;2|^^)(e^|, (93) 

where W 2 is a function of U 2 , and V. 


Exercise 20. Derive Eq. (93) and find the explicit expression of W 2 - 

(3) If —1 is a nondegenerate eigenvalue of both U and V, that is ui = = — 1 and 

U 2 tV 2 ^ ~1, then there are two possibilities: 

(a) If the eigenvectors of U and V belonging to —1 are parallel, that is U com¬ 
mutes with V, then D{tw) = D{tu) = D{tv)- Thus, the only constraint is 
(?l =0 W has the previous form (93). 


Exercise 21. Find how W 2 particularizes in this case to a function of U 2 
and V 2 - 


(b) If the eigenvectors ^ oi U and r] of V belonging to —1 are not parallel, then 
they span the whole space. The constraints (^|'I') = 0 and (r?|'I') = 0 imply 
Dirichlet’s boundary conditions 'h = 0, so that D(tw) = and 


W = U -kV = -I. 


(94) 


(4) Finally, in the case U = —I (or V = —I) then D{tw) = so that 'h = 0 

and 


W = {-I)kU = Uk{-I) = -I. (95) 

Summing up, the emerging boundary conditions are given by the composition 
law (891, where the product * is given by the Cayley transform (92) for U and V 
with no eigenvalues —1 (i.e. free ends 4'). On the other hand, eigenspaces with eigen¬ 
values — 1 are absorbing for the product, that is all constraints on the wave-function 
boundary values 4^ are inherited by W. In particular the Dirichlet boundary con¬ 
ditions —I play the role of an attractor (95), and when U and V have independent 
constraints on 4^, then their composition (94) is Dirichlet, 4> = 0. 


5. Conclusions 

In these lectures we have seen that sometimes examples that are apparently very 
simple may conceal an unexpected rich structure. And in fact they can help us to 
build and understand general schemes. Our goal was to communicate the mathemat¬ 
ical and the physical ideas in the most simple setting, without all the unnecessary 
technical complications that a more realistic model inevitably has. We hope to have 
hit our goal, at least partially. Moreover, we hope that the references we have pro¬ 
vided will stimulate the interest of the reader to go beyond this notes and to get 
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directly in touch with a so active field. For instance, the generalization to higher 
dimensions is nontrivial and introduce new savory ingredients, because in such a 
case the Hilbert space of the boundary is infinite dimensional. But that’s another 
story. 
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